:arXiv:hep-ph/0003136v2 I Jun 2000: 



Nuclear Shadowing and Coherence Length 
for Longitudinal and Transverse Photons 



Boris Kopeliovich 1 ' 3 , Jorg Raufeisen 1 ' 2 , and Alexander Tarasov 1 ' 3 

February 7, 2008 

1 Max-Planck Institut fur Kernphysik, Postfach 103980, 69029 Heidelberg 
2 Institut fur Theor. Physik der Universitat , Philosophenweg 19, 69120 Heidelberg 
3 Joint Institute for Nuclear Research, Dubna, 141980 Moscow Region, Russia 



Abstract 

Motivated by the recent results for DIS off nuclei from the HERMES experiment we 
have performed a systematic study of shadowing for transverse and longitudinal pho- 
tons. We found that the coherence length which controls the onset of nuclear shadow- 
ing at small XBj is much longer for longitudinal than transverse photons, and is much 
shorter for shadowing of gluons. The light-cone Green function approach we apply 
properly treats shadowing in the transition region xbj > 0.01. It also incorporates the 
nonperturbative effects and is legitimate at small Q 2 . We calculate nuclear shadowing 
and compare with data from the HERMES and NMC experiments. Although we ex- 
pect different nuclear shadowing for longitudinal and transverse photons, numerically 
it cannot explain the strong effect observed by the HERMES collaboration. 
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1 Introduction 



Recently the collaboration HERMES released data for shadowing in inclusive positron 
scattering off nuclei at medium high energies and Q 2 . The results expose few unusual fea- 
tures. The cross sections on nuclear targets, nitrogen and helium-3, at small Xsj ~ 0.02 and 
Q 2 < 1 GeV 2 was found to be substantially more shadowed that one could expect extrapo- 
lating available data at higher Q 2 and energies. Unusual is also the observed enhancement 
of shadowing with Q 2 . Interpreted in terms of different shadowing for transverse and lon- 
gitudinal photons it was concluded in [[[]] that a L is enhanced, while <jt is suppressed on 
nitrogen by at least factor of two compared to deuteron target. 

These data drew attention to the fact that very few data are available in this kinematical 
region. Moreover, no reliable theoretical calculations are done yet. The approach based on 
the nonlinear evolution equations ||, ||[ needs knowledge of the nuclear parton distribution 
at a medium-hard scale which is to be guessed, and is anyway outside the kinematical range 
we are interested in. A more promising is the intuitive approach treating nuclear effects in 
the spirit of vector dominance model (VDM) j| as shadowing for the total cross section 
of hadronic fluctuations of the virtual photon (see e.g. in ||). However, the perturbative 
QCD treatment of the photon fluctuation can be applied only at high Q 2 , while VDM is 
sensible only at small Q 2 — > 0. 

A progress was done recently || on extension of perturbative QCD methods to the 
region of small Q 2 where the quarks in photon fluctuations cannot be treated as free. The 
nonperturbative inter-quark interaction was explicitly introduced and new light-cone qq 
wave functions were derived which recover the well know perturbative ones at large Q 2 . 

Nuclear shadowing is controlled by the interplay between two fundamental quantities. 

• The lifetime of photon fluctuations, or coherence time. Namely, shadowing is possi- 
ble only if the coherence time exceeds the mean internucleon spacing in nuclei, and 
shadowing saturates (for a given Fock component) if the coherence time substantially 
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exceeds the nuclear radius. 



Equally important for shadowing is the transverse separation of the qq. In order to be 
shadowed the gg-fluctuation of the photon has to interact with a large cross section. 
As a result of colour transparency JF], |8|, Q , small size dipoles interact only weakly and 
are therefore less shadowed. The dominant contribution to shadowing comes from the 



large aligned jet configurations [10, H of the pair. 



The mean lifetime of a qq fluctuation in vacuum calculated in section 24. turns out to be 
zero for transverse photons. This strange result is a consequence of an incorrect definition. 

In sections [2.2| and |2]^ we propose a more sophisticated treatment of the coherence length 
or the fluctuation lifetime relevant for shadowing. The mean coherence time for the qq Fock 
state is evaluated using the perturbative and nonperturbative wave functions. The salient 
observation is that the coherence length is nearly three times longer for longitudinal than 
for transverse photons. At the same time, both are substantially different from the usual 
prescription l c = {2m^XBj)^ 1 ■ The coherence length is found to vary steeply with Q 2 at 
fixed XBj and small Q 2 . 

The coherence time for a \qqG) Fock component controlling nuclear shadowing for gluons 
is calculated in section |2]4]. It turns out to be much shorter than for \qq) components, 
therefore, onset of gluon shadowing is expected at smaller xsj than for quarks. 

The transition region between no-shadowing at XBj ~ 0.1 and saturated (for the \qq) 
component) shadowing at very small XBj is most difficult for theory. The impact parameter 
representation assigns definite cross sections to the fluctuations, but no definite mass which 
one needs to calculate the phase shift. On the other hand, the eigenstates of the mass matrix 
cannot be associated with any definite cross section. This controversy was settled within 



the light-cone Green function approach JO], [12], I3fl . In section |3| we rely on this approach 
to calculate nuclear shadowing in the kinematical region of the HERMES experiment. The 
nonperturbative light-cone wave functions and the realistic phenomenological dipole cross 
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section are important at low Q 2 . We are unable to reproduce the data from the HERMES 
experiment, although our parameter- free calculations are in a good agreement with NMC 
data. We do not expect any dramatic enhancement of the longitudinal cross section com- 
pared to the transverse one. At the same time, the calculated transverse cross section for 
nitrogen is shadowed by about 20% only, much less than the HERMES data need. 



2 The mean coherence length 

2.1 The lifetime for a perturbative qq fluctuation in vacuum 

A photon of virtuality Q 2 and energy v can develop a hadronic fluctuation for a lifetime, 

2 v P 

1 = = = pimax f i \ 

Lc Q 2 + M 2 x Bj m N Lc ' [) 

where Bjorken xsj = Q 2 /2m a^, M is the effective mass of the fluctuation, factor P = 
(1 + M 2 1 Q 2 )^ 1 , and l™ ax = Xjm^XBj- The usual approximation is to assume that M 2 « Q 2 
since Q 2 is the only large dimensional scale available. In this case P = 1/2. 

The effective mass of a noninteracting q and q is well defined, M 2 = (m 2 + fcfi) / a(l — a), 
where m q and kx and a are the mass, transverse momentum and fraction of the light-cone 
momentum of the photon carried by the quark. Therefore, P has a simple form, 



Q 2 a{l-a) 

P(k T ,a)= fc2+g2 , (2) 



where 

e 2 = a(l - a)Q 2 + m\ . (3) 

To find the mean value of the fluctuation lifetime in vacuum one should average (0) over 
and a weighted with the wave function squared of the fluctuation, 





P(k T ,a) 


* qq / 


\*qq 


* qq / 



(P)vac / T -v* T -v*\ ' (^) 
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The perturbative distribution function for the qq has the form [14, 15, 16 



x6 T ' L X K (er T ) 



2tt 



(5) 



Here \ and x are the spinors of the quark and antiquark respectively. K^erx) is the modified 
Bessel function. The operators T,L have the form, 



T = m q a ■ e + i{\ - 2a) (a ■ n) (e • V T ) + {3 x e) • V T 



O 1 



2 Q a(l — oi) a ■ n , 



(6) 



(7) 



where the two-dimensional operator Vt acts on the transverse coordinate ?t\ n = p/p is a 
unit vector parallel to the photon momentum; e is the polarization vector of the photon. 

The normalization integral in the denominator in the r.h.s. of (^) diverges at large &t 
for transversely polarized photons, therefore we arrive at the unexpected result (P T ) vac = 0. 



2.2 Coherence length in nuclear medium 

This puzzling conclusion can be interpreted as a result overwhelming the fluctuations of a 
transverse photon by heavy qq pairs with very large kx- Such heavy fluctuations indeed have 
a very short lifetime. However, they also have a vanishing transverse size tt ~ 1/&t and 
interaction cross section. Therefore, such fluctuation cannot be resolved by the interaction 
and do not contribute to the DIS cross section. To get a sensible result one should properly 
define the averaging procedure. We are interested in the fluctuations which contribute to 
nuclear shadowing, i.e. they have to interact at least twice. Correspondingly, the averaging 
procedure has to be redefined as, 





P{k T ,a) 


f(l* -> qq)) 


(/(t* qq) 


fd* -» qq)) 



(8) 



where /(7* — > qq) is the amplitude of diffractive dissociation of the virtual photon on a 
nucleon 7* N — > qq N. 
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Thus, one should include in the weight the interaction cross section squared a^ q (r T ,s), 
where is the transverse separation and s = Im^v — Q 2 + m 2 N . Then, the mean value of 
factor P(a, hA reads, 



P 
with 



Jq da J d 2 rid 2 r 2 




Vqq (r 2 , S) P {f 2 ~ fx, O) tyqf (fi, a) (T^ (n, s) 


Jo da J d 2 r 




2 



~ r d 2 k -> 

P (f 2 - fi, a) = J — exp (-i k T ■ (f 2 - f^P (a, k T ). 



Using expression (pj) one obtains for a non interacting gg-pair, 

Q 2 a(l -a) 



P\T2 - f^a) 



2tt 



-K (e\f 2 — r 1 1 ) . 



(9) 



(10) 



ir 



As a simple estimate for the mean value @ one can use the small-r^ approximation for 
the dipole cross section (jq q {r T , s) = C(s) r 2 -,. The factor C(s) does not enter the result since 
it cancels in (^). We obtain for transverse and longitudinal photons, 

)2 Jda(l - a) a ( a 2 + (1 - a) 2 / e 6 + | m 2 q (l - a) a / e 8 ) 



P 



T 



2Q 2 i 



}da{[a 2 + (1 -a) 2 / 



(12) 



P 1 



7 rft I da(l — a) 3 a 3 / e 8 
' v o ' 

8 ]rfa(l-a) 2 a 2 / e6 



(13) 



respectively. 



We calculated the factor (P > ) as function of Q at XBj = 0.01 which is close to the 
minimal value in the HERMES data. We used the effective quark mass m q = 0.2 GeV which 
corresponds to the confinement radius and allows to reproduce data on nuclear shadowing 
[Til . Our results depicted in Fig. |l| by dotted lines are quite different from the naive estimate 
P T ' L = 1/2. Besides, P L turns out to be substantially longer than P T . This indicates that a 
longitudinally polarized photon develops lighter fluctuations than a transverse one. Indeed, 
the effective mass M is maximal for asymmetric pairs, i.e. when a or 1 — a are small. 
However, such fluctuations are suppressed in longitudinal photons by the wave function (0). 



The dependence of I^P T < L j on xsj depicted in Fig. ^| for Q 2 = 4 and AOGeV is rather 
smooth. Therefore, the coherence length varies approximately as l c oc 1/xBj- 

The simple approximation aq q oc r\ is not realistic since nonperturbative effects affect 




0.01 0.1 1 10 100 

Q 2 / GeV 2 



Figure 1: Q 2 dependence of the factor (P) = l c jl™ ax defined in (fj) at 
XBj = 0.01. for qq fluctuations of transverse and longitudinal photons, 
and for qqG fluctuation, from the top to bottom, respectively. Dotted 
curves correspond to calculations with perturbative wave functions and 
an approximate dipole cross section oc r\. Dashed curves are the same, 
except the realistic parameterization ([T^j . The solid curves show the 



most realistic case based on the nonperturvative wave function (|25). The 



coherence length for gluons calculated in sec. is also shown. 
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Figure 2: XBj dependence of the factor (P T < L j and \P G J defined in 
(W corresponding to the coherence length for shadowing of transverse 
and longitudinal photons and gluon shadowing, respectively. Solid and 
dashed curves correspond to Q 2 = 4 and 40 GeV 2 . 



the large-7Y behaviour. Motivated by the phenomenon of confinement one should expect 
that gluons cannot propagate far away and the cross section should level off at large tt- 
We use the modification H of the energy dependent phenomenological dipole cross section 



cr^{r Tl s) suggested in [17 



^-0t,s) = ob(s) 



qq 



1 — exp 



rg(s) 



(14) 



where r (s) = 0.88 (s/sq) 0,14 fm, Sq = 1000 GeV 2 . This cross section is proportional to r\ 
at small ~~ * 0, but is constant at large Tt- The energy dependence correlates with r^, at 
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small 7Y the dipole cross section rises steeper with energy than at large separations: 

3r 2 (s)\ 



crZ(s) 1 + 



( r ch)TT, 



ch 



0.44 fm 2 



(15) 



where 



aZ(s) = 23.6 {a I s ) om mb. (16) 

With this choice of cr^(s) one automatically reproduces the total cross section for pion 
proton scattering, while the parameterization from |T7[ cannot be applied to hadronic cross 



sections. Thus, cross section is better designed for low and medium large Q < 10 — 
20GeV 2 , while at high Q 2 the parameterization |T7| works better. 
Eq. (0) can be represented in the form, 



P 



T,L 



N 



T,L 



(17) 



D T,L ' 

The angular integrations in (1.7) for the denominators D T ' L are trivial and for the numer 
ators N T,L one uses the relation [|l~8|j , 



K (e |n - r 2 |) = K Q (er>) J (er<) + 2 £ e 4m ^ m (er>) J m (er 



m=l 



where r> = max (7*1,^), r< = min(ri,r2), cos0 = n • ^2/(^1^2) and I m {z) are the modified 
Bessel functions of first kind (Bessel function of imaginary variable). It is clear from this 
relation that after angular integration only one term in the sum gives a non-vanishing 
contribution. We finally obtain for transverse photons, 



Ni 



+ 



2 Q 2 J data (1 - a) J dr 2 r 2 J dr x r 1 jm 2 K% (er 2 ) K (en) I (eri) 
00 

e 2 K\ (er 2 ) K x {er,) h (en)} o» (s, n) < (s, r 2 ) , (19) 



a 2 + (l-af 



D T V = j da j dr r {m 2 , K 2 (er) + 


and for longitudinal photons, 



a 2 + (1 -a) 2 



Kl(er)} 



(20) 



N$ = 2 Q 2 J da a 3 (1 - af J dr 2 r 2 J dr x r x K 2 (er 2 ) K (en) I (en) (s, n) < (s, r 2 ) 



(21) 
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Dp = J da J drra 2 (1 - a) 2 K 2 (sr) (s, rf . (22) 
o 

The factor (p T ' L {x ) Q 2 )^j calculated in this way is depicted by dashed lines in Fig. [I] as 
function of Q 2 at XBj = 0.01. It is not much different from the previous simplified estimate 
demonstrating low sensitivity to the form of the dipole cross section. 

It is instructive to compare our predictions with the VDM which is usually supposed to 
dominate at small Q 2 < m 2 . The corresponding coherence length l^ DM is given by ([!]) with 
M = m p . The ratio of Ij calculated with the nonperturbative wave function to l^ DM as 
function of Q 2 is shown by solid curve in Fig. |3]. It demonstrates an unexpectedly precocious 
violation of VDM at quite low Q 2 . We also calculated with the perturbative wave function, 
but with a massive quark. With m q = 200 MeV it mimics the nonperturbative effects quite 
well, as one can see from Fig. ^. 




0.4 - 
0.2 - 
0.0 - 



1 2 

Q 2 / GeV 2 

Figure 3: Q 2 dependence of ratio of (/J) calculated with Eq. (||) and 
m q = 200 MeV to l v c DM calculated with Eq. (fj) and M = m p . 
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2.3 Propagation of interacting qq 

Although the quarks should be treated perturbatively as nearly massless, at the endpoints 
a or 1 — a —>■ the mean qq transverse separation ~ 1/e becomes huge ~ l/m q . 
This contradicts the concept of confinement and should be regularized either by explicit 
introduction of a nonperturbative interaction between q and q f|, or as is done below, 
introducing an effective quark mass. 

Apparently, it is not legitimate to use the perturbative qq wave functions (|5|) - (Q) at 
low Q 2 where the nonperturbative interaction between q and q becomes important. This 
interaction squeezes the qq wave packet, i.e. increases the intrinsic transverse momentum 
and the effective mass of the pair. By contributing to the effective mass of the qq, the 
nonperturbative interaction breaks down the validity of the kinematical expression for M or 
Eq. (|2]). Even at high Q 2 at the endpoints a or 1 — a — > the mean qq transverse separation 
r T ~ 1/e becomes huge ~ 1/ m q and nonperturbative corrections may be important. One can 
try to mimic these effects by an effective quark mass, as is done above, but one never knows 
how good this recipe is. To take the nonperturbative effects into account we use the light- 



cone Green function formalism suggested in || generalizing the perturbative description [|TT 
of nuclear shadowing in DIS. 

Propagation of an interacting qq pair in vacuum with initial separation f% at the point 
with longitudinal coordinate Z\ up to point z 2 with final separation r 2 is described by a Green 
function Gq q (r 2 , z 2 \r\,z\) which is a solution of a two dimensional Schrodinger equation, 



dz 2 2fo(l-a) 



The nonperturbative oscillator potential contains function a(a) = ao + a\ a (1 — a) with pa- 
rameters ao and a± fixed by data ||. Data on total photoabsorption cross section, diffraction 
and nuclear shadowing are well described with, 

a 2 (a) = v 1A \ll2MeV) 2 + (1 - v) 1A5 (IQhMeVfa (1 - a) , (24) 

where v is any number satisfying < v < 1. 
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The Green function allows to calculate the nonperturbative wave function for a qq fluc- 
tuation, 



\l/-' L ff a) -- 
qq 1 ' ' 4vrpa(l -a) 



iZ q y/oT e 



dz 1 (x0 1 ' L X )G™ c (z 1} r 1 ;z 2} r 2 



ri=0; T2=f 



(25) 



In the limit of vanishing interaction, a(a) — > the wave functions (j25l) recovers the pertur- 
bative ones Eq. (H). 

Note that in kr representation the free Green function G^ (zi, fi; z 2 ,f 2 ) integrated over 
longitudinal coordinate is simply related to the coherence length (JTJ), if one performs an 
analytic continuation to imaginary time, z — > —iz 



dz 2 G^z^fi, z 2 ,f 2 ) 



d 2 ki 



exp 



-i k T ■ (r 2 - ri) l c (k T , a) ■ 



(26) 



(2tt) 2 

This is easily generalized to include the nonperturbative interaction. Then, making use of 
this relation one can switch in (|S|) to vt representation and express the mean coherence 
length via the Green function. Then, we arrive at new expressions for the functions N T,L 
and D T ' L in (IT! 



N T,L 



da I d\ x d z r 2 ^ q L (f 2 , a) 



dz 2 G v q a g c {f u Zl ; f 2 , z 2 ) 



D 



T,L 



da / d 2 r 



^ L (ra)a»(r,s) 



where the nonperturbative qq wave functions are given by fl25|) . 

For a harmonic oscillator potential the Green function is known analytically, 



(f 2 , Zl - r^z,) 



a 2 (a) 



27rsinh (u Az) 



exp 



e 2 Az 



2 v a(l — a) 



x exp 



a 2 (a) 



r i + r 2 ) coth (uj Az) 



2fi • f 2 
sinh (u Az) 



where Az = z 2 — z\ and 
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a(a) 2 
vail — a) 



(27) 
(28) 



(29) 



(30) 
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is the oscillator frequency. 

The light-cone wave functions fl25f ) modified by the interaction can now be calculated 
explicitly [§, 



g9 V~) "j 



Z 2 a 



m 2 $* (e, A, $ (e, A, r 2 ) 



1 _ 2a (1 - a ) $* ( e , A, fx) • $ x (e, A, f 2 ) 



^ (e, A, f 2 ) = Z^a< em 4QV (1 - a)' $q (e, A, ri) $ (e, A, f 2 ) , 



where the parameter 



A 



2a 2 (a) 



-2 ' 



describes the strength of the interaction; and the functions $o,i read, 



$* (e,A,fi)$ (e,X,f 2 ) = -—2 

(4tt) 



c/m (it 



A 2 



sinh (At) sinh (Am) 



x exp 



coth (At) - t - coth (Am) - u 



(31) 
(32) 

(33) 



(34) 



(e,A,ri) (e,A,r 2 ) 



1 n • r 2 

(2vr) 2 r?r| 



G?W (it 



x exp 



' 1 coth (At) - t - bill co th (Am) - u 



(35) 

4 4 
It is easy to verify that in the limit of vanishing interaction, A — > 0, the nonperturbative 
wave-functions reduce to the perturbative ones. Compared to the expression for $i in || , we 
have integrated by parts over the parameter u, or t respectively. This considerably simplifies 
the expression. 

Now we have all ingredients which are necessary to calculate flT7|). Two from the eight 
remaining integrations, over the angles, can be performed analytically. We obtain, 



N 1 



m^XBj j dot I driridr 2 r 2 



o 

N I \ N / \ ^ v'-V 



dAz 



x Ii 



qq v t, - qq 

a 2 (a) r x r 2 



sinh (uo Az) 



exp 



^(e,A,r 2 ) fygq (e, A, ri) 



2 i/ ct(l — a) 



sinh (a; A2) 



exp 



a 2 (a) 



(r\ + r 2 ) coth (00 Az) 



(36) 



14 



N 1 



m N x B j 



da 



x <J Q g(r 2 ,s) a,Jr u s) 



dr± r± dr 2 r 2 
a 2 (a) 



sinh (lu Az) 



dAz 
exp 



e 2 Az 
2 v a(l — a) 



*L A, ri) 



a 2 (a) r\ r 2 
sinh (u Az ) 

D L,T = 



exp 



da 



a 2 (a) 



r\ + r 2 J coth (a; Az) 







tt^ (e, A, f) < (r, S ) 



</</ 



(37) 



(38) 



For a dipole cross section that levels off like (|14D at large separations the integrations over 
7*1 and r 2 can also be done analytically. However, we prefer to work with the more general 
expressions that hold for arbitrary cr^ (s, r), as long as it depends only on the modulus of r. 
We perform the remaining integrations numerically. The results for l^ ,L (x,Q 2 ) are shown 
by solid curves in Figs. [1| and 0. 



2.4 Coherence length for gluon shadowing 

Shadowing in the nuclear gluon distributing function at small XBj which looks like gluon 
fusion GG — » G in the infinite momentum frame of the nucleus, should be treated in the 
rest frame of the nucleus as shadowing for the Fock components of the photon containing 
gluons. Indeed, the first shadowing term contains double scattering of the projectile gluon 
via exchange of two t-channel gluons, which is the same Feynman graph as gluon fusion. 
Besides, both correspond to the triple-Pomeron term in diffraction which controls shadowing. 

The lowest gluonic Fock component is the \qqG). The coherence length relevant to 
shadowing depends according to ([[]) on the effective mass of the \qqG) which should be 
expected to be heavier than that for a \qq), and even more for higher Fock components. 
Correspondingly, the coherence length (if) should be shorter and a onset of gluon shadowing 
is expected to start at smaller XBj- 

For this coherence length one can use the same Eq. flU), but with the effective mass, 

P M 2 

Ml G = n kT r + T^, (39) 
a G (l - a G ) 1 - a G 
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where «g is the fraction of the photon momentum carried by the gluon, and M qq is the 
effective mass of the qq pair. This formula is, however, valid only in the perturbative limit. 
It is apparently affected by the nonperturbative interaction of gluons which was found in |j| 
to be much stronger than that for a qq. Since this interaction may substantially modify the 
effective mass M qqG we switch to the formalism of Green function described above, which 
recovers Eq. (|39|) in the limit of high Q 2 . 

We treat gluons as massless and transverse. For factor P defined in ([!]) one can write, 

= w • < 40 > 

where 

N G = m N x B j J d 2 r 1G d 2 r lqq d 2 r 2G d 2 r 2qq da q dln(a G ) ^\ qG {f 2G ,r 2gq , a q ,a G ) 

x \fdz 2 G qqG {r 2G ,r 2qq ,z 2 \ri G ,r lqq ,z{^ ^q q G(ri G ,r lqq ,a q ,a G ) (41) 

D° = J <er lc *r» <fr 20 d % 5 la, i ln(a o) §*, G (f 2G , f m , a„a G ) 

x 5 {2) (r 2G - rio) S {2) (f 2qq - f lqq ) ^ qqG (fio, r lqq , a q , a G ) (42) 

Here we have introduced the Jacobi variables, f qq = R q — R q and r G = R G — (a q R q + 
a q R q )/(a q + a q ). R G , q , q are the position vectors of the gluon, the quark and the antiquark 
in the transverse plane and otG,q,q are the longitudinal momentum fractions. 

The Green function describing propagation of the qqG system satisfies the time evolution 
equation 0, 

' d Q 2 a q + a q A ± (rag) " 

a 7T - + o A -L l r qq + 7, n \ ~ V \ r 2G, r 2qq, ®q, «G, %2 

oz 2 2v 2va q a q 2va G (1 — a G ) 

x G qqG (r 2G , f 2qq , z 2 ; f 1G , r lqq , z x ) = 8{z 2 -z 1 ) 5 {2) (f 2G - f XG ) 5 {2) (f 2qq - r lqq ) .(43) 

In order to calculate the coherence length relevant to shadowing, we employ the ampli- 
tude for diffractive dissociation 7* — > qqG, which is the qqG wave function weighted by the 
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cross section 



^qqG (r G , fgg, Ot q , <X G ) = tfj^ (f qq , a, 



Oin ~\~ OLq 



On 



OLq ~\~ Oin 



gg 



9 




8 


°gg\ 



a,-, 



gg 



gg 



OCr, 



gg 



r gq ) a qq i x i r gg) 



(44) 



\ C^q Oiq 

As different from the case of the \qq) Fock state, where perturbative QCD can be safely 
applied at high Q 2 , the nonperturbative effects remain important for the \qqG) component 
even for highly virtual photons. High Q 2 squeezes the qq pair down to a size ~ 1/Q, however 
the mean quark-gluon separation at a G <C 1 depends on the strength of gluon interaction 
which is characterized in this limit by the parameter 6 ~ 0.65 GeV |]]. For Q 2 3> 6q the qq 
is small, r| <C r G , and one can treat the qqG system as a color octet-octet dipole, i.e., 

G q q G (r 2G ,r2qq,Z2;r 1G ,r lq q,z 1 ) =$> G q q(r 2 qq,z 2 ;r lq q,z 1 ) G GG (f 2G ,z 2 ;f 1G ,z 1 ) . (45) 

Such a Green function G GG satisfies the simple evolution equation ||, 

G GG (f 2 G,Z2;r 1G ,Zi) 



d Q 2 + A ± (r 2G ) 



K r 2G 



dz 2 2v 2ua G (1 — a G ) 2ua G (1 — a G ) 

= 5{z 2 - z{) S {2) {r 2G - t X g) 



(46) 



Correspondingly, the modified qqG wave function simplifies too, 



^qqG (r G , f qq , U q , Oq) => ~ * g (f qq , O q ) f G ■ V * qG (f G ) O GG (x, T G ) 



(47) 



where the nonperturbative quark-gluon wave function has a form 



T / - \ ! • t / \ 2z [aj e ■ r G 
^qG (r G ) = km ty qG (a G , r G ) = -— x j — — exp I —jr r G . , 



b 2 



7T V 3 r 



G 



and the color-octet dipole cross section reads, 



(4* 



°gg K x , r G) = t cr qq - [x, r G ) . 



(49) 
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With the approximations given above, the factor (P G ) in ([!]) for the gluon coherence 
length is calculated in Appendix A and has the form, 

^vnax 



G 



P 



where 



3 ln(a% ax /a% in ) 5 



dS 



s 



8(1 + 5) 8(1 + 35) 5 2 - 1 



(50) 



d\nT(x) 2% 

^ ) = ^*— ' 6 = o^- (51) 

Both the numerator and denominator in (|50"D diverge logarithmically for a™" — > 0, as 
it is characteristic for radiation of vector bosons. To find an appropriate lower cut off, note 
that the mass of the qqG system is approximately given by 

< G -^ + Q 2 , (52) 

where we used eq. (|39|) with (kj) ~ fr 2 ,. We demand that M 2 ^ G < 0.2s which leads to 
a™ m = 26g/ (0.2s — Q 2 ). Furthermore we work in the approximation of olq <C 1 and we also 
have to choose an upper cut off. We use 

2bl 2hl . . 

< «g < 7^ (53) 



0.2s - q 2 - Q 2 

which means that we take only masses 2Q 2 < M 2 ^ G < 0.2 s into account. The two limits 
become equal at %Bj ~ 0.1. 

Our results for (P G ) = (l G )/l™ ax are depicted in Fig. |I[ With approximations made 
above we cannot cover the low Q 2 region and perform calculations at Q 2 > lGeV 2 . The 
found coherence length is much shorter than both Ij and if for \qq) fluctuations. This 
conclusion corresponds to delayed onset of gluon shadowing shifted to smaller xsj predicted 
in!- 



3 Shadowing for longitudinal and transverse photons 
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3.1 fr L /cr T on a nucleon target 

As soon as realistic wave functions for qq fluctuations including the nonperturbative effects 

are available, as well as the energy dependent phenomenological dipole cross section, we are 

in position to calculate the longitudinal and transverse cross sections for a proton target 

covering also the region of small Q 2 , 

1 

a T,L = J da J d 2 r 
o 

The results of calculations for the ratio <Jl/&t is shown in Fig. |] by solid curve as function 
of Q 2 at XBj = 0.01 which is about the lowest value of XBj in the HERMES data. As one 
could expect the ration vanishes towards Q 2 = 0, however is nearly constant down to very 
small Q 2 ~ 0.3 GeV 2 . To see how well an effective quark mass can mimic the effect of the 
nonperturbative interactions we have also performed calculations with the perturbative wave 
functions and m q = 200 MeV and plotted the results by dashed curve in Fig. [|. Comparison 
demonstrates a substantial difference at small Q 2 . 

3.2 Nuclear targets 

Although the coherence length is an important characteristics for shadowing, it is not suf- 
ficient to predict nuclear effects in the structure function. Shadowing for parton densities 
at small XBj in the nuclear structure function which is defined in the infinite momentum 
frame originates from the nonlinear effect of parton fusion in the evolution equation |2|, ||. 
Although the partonic treatment varies from to another reference frame, all observables 
including shadowing are Lorentz invariant. In the rest frame of the nucleus shadowing in 
the total virtual photoabsorption cross section a] ot A (or the structure function F^) can be 
decomposed over different Fock components of the photon, 

°tt = A(T tot N ~ Aa tot (qq) - Aa tot (qqG) - Aa tot (qq2G) - ... (55) 

According to above calculations and Fig. [| the coherence length corresponding to gluon 
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T,L 
99 



fa) 



a 



99 



r,s) 



(54) 



0.4 



0.3 



D 
D 



0.2 



0.1 



0.0 




Q 2 / GeV 2 



Figure 4: Ratio of longitudinal to transverse photoabsorption cross sec- 
tions as function of Q 2 at %Bj = 0.01. The solid curve is calculated 



with (|5i|) and the wave functions including the nonperturbative ef- 
fects, while for the dashed curve the perturbative wave functions with 
m q = 200 MeV are used. The NMC data points JT^/ correspond to 
x Bj = 0.008, 0.0125, 0.0175 from smaller to higher Q 2 . 



shadowing is rather small compared to the mean internucleon spacing in a nucleus in the 
kinematical region we are interested in. Therefore, we hold only the first two terms in the 
r.h.s. of Eq. ( j55|) , which can be represented for transverse and longitudinal photons as, 

i 



(vlof) ' = A i^tot 1 ) ' ~ 7^R e J d 2 b J da j dz x j dz 2 J d\\ J d 2 r 2 



(56) 



Z\ 



^ T qq L (e, \ r 2)\ Pa {b, z 2 ) ( s > r 2) G (f 2 , z 2 I fi, Zl ) p A {b, z x ) (s, n) ^ q L (e, A, n) 
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where z) is the nuclear density dependent on impact parameter b and longitudinal 
coordinate z. The nonperturbative wave functions for the qq component of the photon 
are defined in Eqs. fl3T|)-(|32"D. The Green function G (r 2 , z 2 \ fx, Z\) describes propagation of 
a nonperturbatively interacting qq pair in an absorptive medium. It fulfils the evolution 
equation, 



. d A ± (r 2 ) - e 2 i , N . . a 4 (a) r\ 

^ + 77^7? T + «Pa (b, z 2 ) a" (s, r 2 ) - 



G(f 2 ,z 2 \ri,zi) 



dz 2 2ua (1- a) 2 r " v ' A ' qq v ' zy 2ua (1 - a) 

t5(z 2 -z 1 )5 {2) (f 2 -f 1 ) , (57) 



where a(a) and A were introduced in Eqs. (^1) and ( P5[ ) respectively. The third term in 
the l.h.s. of (|5TD describes absorption of the gg pair in the medium of density pa (b, z) with 
cross section a^ q (s,r). 

At small XBj when the coherence length substantially exceeds, l^ ,L ^> (z 2 — z\) (the 
nuclear radius) the solution of Eq. ( p^) very much simplifies, G(r 2 , z 2 | r*i, zi) oc [r 2 - ri ), 
i.e. Lorentz time dilation "freezes" variation of transverse qq separation. 

Correspondingly, the total cross section gets a simple form M, |16|], 



KilL = *jdajd 2 r \ti%> (er)\ 2 J d 2 b 



exp ^AZT (6) 



(5? 



T(6) = / dzp A (b,z) (59) 



where 



is the thickness function of the nucleus. 

Eq. (|57D has an explicit analytical solution only if the dipole cross section cr^(r) = C r 2 
and the nuclear density is constant pA{b,z) = po- Such an approximation has a reasonable 
accuracy, especially for heavy nuclei. Nevertheless, it can be even more precise if one makes 
use of the fact that the asymptotic expression (^Sj) is easily calculated with exact (realistic) 
dipole cross section and nuclear density. We need to use the solution of Eq. fl5"T|) only in the 
transition region from no-shadowing to a fully developed shadowing given by (0). First 
of all, we fixed factor C in the simplified the dipole cross section demanding to have the 
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same asymptotic shadowing (|58[) as with realistic one given by (14). This was done with 
the realistic Woods-Saxon form for nuclear density and separately for transverse and 
longitudinal photons and for each value of a. Then we switched to a constant nuclear density 
Po, demanding it to lead to the same asymptotic shadowing in (|5^) as with the realistic one. 
We have checked that the found value of p is practically independent of the value of the 
cross section in the interval 1 — 50 mb. 

First of all we have checked our formalism comparing with the NMC results for shad- 
owing in the nuclear structure function. Fig. |5] demonstrates the data [21] for tin to carbon 
ratio of proton structure functions depicted by full circles. As was pointed out above gluon 
shadowing is negligibly small at XBj > 0.01 which covers the whole range on the NMC exper- 
iment. We performed calculations with parameter v = 0.5 in fl24j), but we have checked that 
the results are independent of v. Although the calculations are parameter free, agreement 
is pretty good. 

Eventually, we are able to provide predictions for the kinematical range of HERMES. The 
ratio of the virtual photoabsorption cross sections for nitrogen to hydrogen at XBj = 0.01 
versus Q 2 is plotted in Fig. |. The solid and dashed curves correspond to the nonperturba- 
tive and perturbative wave functions respectively. A salient feature of these predictions is 
shadowing vanishing towards Q 2 = 0. This observation does not contradict the well known 
fact of shadowing for real photo absorption, but is a rather simple consequence of kinemat- 
ics. The colliding energy s = Q 2 XBj vanishes along with Q 2 at fixed XBj- This example 
shows that xbj is a poor variable at small Q 2 and should be replaced by s. This kinematical 
effect may be partially responsible for the unusual enhancement of nuclear shadowing with 
rising Q 2 detected in the HERMES experiment |l| . Our calculated nitrogen to proton ratio 
of the cross sections reaches its minimum at Q 2 ~ 0.2 GeV 2 and then smoothly rises with 
Q 2 - 

The most striking feature of the HERMES data for shadowing is a dramatically rising 
R = ctl/gt ratio on nitrogen compared to proton target at Q 2 < 1 GeV 2 [l|. Our predictions 
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Figure 5: Q 2 and XBj dependence of structure function ratio for tin to 
carbon. Full circles show the NMC data The crosses show the 

results of our calculations for the same kinematics. 



for R N /R p are plotted in Fig. [7] versus Q 2 at XBj = 0.01 where the experimental ratio 
Rn / Rp ~ 5. Apparently, we do not expect any remarkable effect. Moreover, this ratio does 
not change much from proton to nitrogen in spite of the much longer coherence length for 
longitudinal photons predicted above. However, this effect leading to a stronger shadowing 
for <7l is compensated by the fact that longitudinal photons develop fluctuations of a smaller 
size compared to that in transverse photons, i. e. they are less shadowed. Note that nuclear 
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0.7 1 1 

1 2 

Q 2 / GeV 2 

Figure 6: The shadowing ratio for nitrogen over proton at low Q 2 and 
x Bj = 0.01. Shadowing disappears as Q 2 — > 0, because of the vanishing 
coherence length. 

effects for R were estimated previously in |22| where many of important ingredients of 
present approach were missed. 

4 Summary 

Basing on the light-cone Green function approach we provide predictions for nuclear shad- 
owing in the most difficult for calculation region of medium-small Xsj > 0.01. Since the non- 
perturbative effects are included we are also predict shadowing down to small Q 2 < 1 GeV 2 
which is the kinematical region where the HERMES experiment discovered unusual shad- 



24 
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0.8 1 1 

1 2 

Q 2 / GeV 2 

Figure 7: Ratio of longitudinal to transverse fractions of the cross section 
R = (T L /a T on nitrogen to proton targets. The solid and dashed curves 
correspond to the photon wave function with and without nonperturba- 
tive effects. 

owing effects. 

We found that the coherence length which controls shadowing is nearly three time longer 
for longitudinal than for transverse photons, and it is very much different from what is 
suggested by the widely accepted approximation. 

Using the Green-function light-cone approach including the nonperturbative effects we 
calculated nuclear shadowing for longitudinal and transverse photons. Although the pre- 
dicted nuclear shadowing exposes interesting effects in the region of small Q 2 , we are unable 
to explain the dramatic phenomena detected in the HERMES experiment. 
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We suppose that our results provide a reliable base line for nuclear effects in this region. 
The dramatic effects revealed by the HERMES experiment probably cannot be explained 
without involving a new nonstandard dynamics (e.g. see in f23|). 
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Appendix A Calculation of the coherence length for 

the qqG-state 



The calculations are somewhat more cumbersome for the gluon coherence length. We con- 
sider the case Q 2 3> With the approximations (fl"5| ) to (47) we obtain for the denominator 



of (H) 



D 



G 



d 2 r G d 2 r q q I da q 



da 



G 



a G 



\1/ ( n CqqOLqq ) 



°qq 0, T C) 



r q q-V^ qG (r G ) . (A.l) 



This integral diverges logarithmically for a G m — > 0. To find an appropriate lower cut off, 
note that the mass of the qqG system is approximately given by 

2b 2 



Kg 



(A.2) 
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We demand that M 2 qG < 0.2s which leads to a™ = 2b 2 / (0.2s — Q 2 ). Furthermore we work 
in the approximation of «g <C 1 and we also have to choose an upper cut off. We use 

2b 2 2b 2 
° < oca < (A.3) 



0.2s -Q 2 ~ " ~ Q 2 

The two limits become equal at XBj ~ 0.1. In our further calculation of D G we do the 

qqi r qqj ~^ T qq^ij 

2 6 ^ 

D ° = (^\/y) 7^ 4 ® 2n J d2rcdr ™ r3qq L daq / ~^ K ^ £rq ^ a ^ 1 ~ a ^ 



replacement r qqi r qq ,-. — > ri?=&j and perform the derivative. This yields 



x 



(2tt)< 

2 _ fe2 . 2 / 2 2/- 



2 



<(s,r G ) e -Vo + o +ft g, (A _ 4) 



4 1 2 



For the integration over r qq we use the integral representation 

r, / X 1 f°° dt ( X 2 \ . . . 

*.(*) = 5 j( jexp(-t- ¥ j (A.5) 
of the modified Bessel function and obtain 

2 



y o dr qq r 3 qq -K 2 (er qq ) = —. (A.6) 



Thus we have for the denominator 



32cwa,, /■«>, r , ,i2 (2 2b 2 



° u - w 111 w h iTGTa w < s > r °)J" «■** U + ^ + "°J • (A7) 

Now we restrict ourselves to the a dipole cross section of the form 

a^(s,r G ) = C(s)r 2 G (A.8) 

and perform the last integration with the result 

^ r 32a em a s C 2 (s) , a^ ax 

Note that the factor C(s) will drop out, when one takes the ratio < P G >= N G /D G . 
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Next we calculate the numerator 
2 [a~7\ 2 6a 



N G = m N x B j \ —\/ — 



vr V 3 ) (2tt) 



.. m a x 

l G 

m AQ 2 J d 2 r 1G d 2 r 2G d 2 r qq J da q J 



da 



a G 



G 2/i \2 

-a q (l-a q ) 



x K (er qq ) a qq (s, r XG ) a qq - (s, r 2G 

x / dAzG GG (f 2G ,f 1G ,Az) 
Jo 



r qq " V riG „ e 2 



r <?<? ' r 2G 2 e 2 



2G 



(A.10) 



where 

Cgg (r 2G ,ri G , Az) 



b\e 



2tc sinh feu Az 



exp ^ - 



h 2 

°0 



(/■ 2 G + r 2 2 G ) cth^Az)--^ 



2fi G • f 2G 



sinh (uuAz) 



(All) 



is the solution of (£46|) with c<j = b\j {va G ). 



Again we can do the replacement r qqi r qq j — > r 2 ^^, perform the derivatives, sum over 
gluon polarizations and use ( |A.6| ) to obtain 



, tG ^a em a s f 

N = mNXBj ^&J 



8a em a s i o ,o , . 

' a r\ G d r 2G aAz 



, max 
G 



da G N , , K , , b 2 e~ ^ 



a G 



^ (S, n G ) Oqq (S, r 2G ) 



sinh (uiAz) 



\ {riG ■ r 2G f h4 (fi G ■ f 2G f nh2 (Fig ■ r 2G ) 2 , 2 (f 1G ■ r 2G f 

4 3 ' h °0 ~2 „2 h ZO ~2 „4 h ZU 0~ 



4 4 
' IG' 2G 

/) 2 /) 2 



' 1G' 2G 



r 2 r 4 
' 1G'2G 



4 2 



1G' 2G 



2 2 2 2 

' IG ' IG '1G'2G 



exp (-/3 (r\ G + r 2 G ) + 27r 1G • r 2G 



where 



b 2 

(3 = -2(i + C oth (uAz)) 
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fo 2 



(A.12) 

(A.13) 
(A.14) 



2 sinh [uAz) 

With a cross section like a qq (s, r) = C(s)r 2 the integrations over r G are easily performed 
with the result 



N G = m N x Bj 



% a em a s C 2 (s) 
3tt 2 Q 2 6 2 



dAz 



max _ n 



10 



12 



sinh 2 (uAz) (1 + cth (a;Az)) c 



a G sinh (cuAz) [ (1 + coth (uAz)Y 
- 8 sinh 2 wAz In * '- 



(A.15) 
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8a Pm a K C 2 (s) r 1 



3iT 2 Q 2 b 2 



J^dy J — {Ay^ 2 (l - y 2 ) In (l - y 2 ) + (l - y 2 ) + I2y^ +2 } . 

(A.16) 



In the last step we have introduced the new variables y = e and 5 = 2b\j (Q 2 Qg)- 
According to (|A.3| ) the limits for 5 are 



l<i<0.2^-l. 



(A.17) 



For the integral containing the logarithm it is convenient to do one more substitution, 
x = y 2 . Then one finds 



f 1 dyy^ 2 (l - y 2 ) In (l - y 2 ) 

J 



(A.18) 



- lim [ dxx^s 2 (l — x) 1+r ' 
2 >?->o dr] Jo v 1 



Now only one integration is left in N G 



iV : 



grnax 

a __ 64a em a s C 2 (s) r d5_ 
3n 2 Q 2 b 2 J 5 



7 5 ( , . , ,/3 1 

[i+8) 1 8(iT^)"^rr (2) "^b + 27 



+ 



(A.20) 



and we end up with the result fl50D for the factor (P G ) = N G /D G at Q 2 > 6q. 
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